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In  an  investigation  of  the  mechanism  responsible  for  the  transport  of  water  through  Si02  thin 
films,  we  have  performed  tracer  diffusion  measurements  involving  network  lsO  demonstrating 
the  importance  of  oxygen  exchange  between  the  Si02  network  and  molecularly  dissolved  water. 
We  have  found  that  in  the  presence  of  water,  bound  network  oxygen  diffuses  through  Si02  as  a 
constituent  of  molecularly  dissolved  water.  Employing  methods  common  to  state-of-the-art 
semiconductor  technology,  the  central  region  within  a  thermal  oxide  layer  grown  on  silicon  was 
enriched  with  immobile  lsO  by  ion  implantation.  After  heating  in  atmospheres  with  different 
water  contents,  the  extent  of  ,80  diffusion  was  determined  by  observing  changes  in  the 
concentration  profile  of  implanted  ,90  by  means  of  nuclear  resonance  profiling  making  use  of  the 
629-keV  resonance  of  the  ,80 (p,a)  15N  reaction.  Diffusions  conducted  in  steam  at  1  atm  (at 
temperatures  as  low  as  250  C)  showed  an  activation  energy  of  about  1 7  keal/mol,  which  is  close  to 
that  originally  measured  for  water  permeation  in  Si02.  Diffusions  in  both  air  and  dry  nitrogen 
showed  a  similar  activation  energy,  with  respective  pre-exponential  factors  two  and  three  orders 
of  magnitude  below  the  steam  value.  Diffusion  in  low-pressure  water  vapor  showed  a  clearly 
linear  dependence  on  gas  phase  water  concentration  down  to  80  ppm.  These  results  are  consistent 
with  a  model  of  water  diffusion  in  Si02  in  which  the  diffusion  mechanism  is  the  interstitial 
transport  of  dissolved  molecular  water  accompanied  by  a  reversible  reaction  with  silicon-oxygen 
bonds  in  tie  network. 

PACS  numbers:  66.30.  —  h,  81.60.Dq 
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I.  INTRODUCTION 

The  formation  of  thin  insulating  films  on  the  surface  of 
semiconductor  materials  is  an  important  part  of  modem 
electronic  technology.  Because  of  its  commercial  impor¬ 
tance,  various  means  of  film  growth  have  been  intensively 
studied,1  especially  for  high-quality  gate  insulators  on  sili¬ 
con  surfaces.  For  that  application  a  more  or  less  standard 
practice  of  thermal  oxidation  in  dry  oxygen  or  steam  atmos¬ 
pheres  has  evolved,  by  which  oxide  films  are  routinely  pro¬ 
duced  satisfying  the  stringent  requirements  of  integrated  cir- 
cuit  devices.  However,  the  ever-increasing  complexity  and 
Q»  operating  speed  of  such  devices  requires  ever  smaller  and 
O  more  closely  packed  component  structures.  The  scale  of 
~  sizes  has  now  reached  a  point  where  basic  changes  are  re- 
I  |  |  quired  in  manufacturing  methods.  For  example,  the  need  is 
l  foreseen  to  hold  temperatures  in  all  steps  of  a  processing 
^  sequence  low  enough  to  prevent  significant  migration  of  do¬ 
pants.  Virtually  all  of  the  many  studies  of  thermal  oxidation 
have  been  device  oriented.  Physical  models  of  the  oxidation 
process  have  received  little  attention,  especially  for  steam 
oxidation;  only  now  are  models  being  tested  with  sufficient 
care  to  generate  a  consensus  of  acceptance.  For  the  future,  a 
greater  understanding  of  the  mechanism  responsible  for  oxi¬ 
dation  in  the  presence  of  water  is  imperative,  especially  in  the 
new  low-temperature  regime,  where  it  now  appears  that  pro¬ 
cessing  of  the  coming  generation  of  devices  will  take  place. 

The  main  object  of  our  work  is  to  clarify  the  mechanism 
of  Si02  film  growth  in  atmospheres  containing  water. 
Through  determination  of  the  kinetics  of  network  oxygen 


exchange,  we  provide  a  direct  experimental  demonstration 
of  the  validity  of  a  recently  proposed  model  of  the  oxidation 
process.2  According  to  this  model,  growth  proceeds  by  in¬ 
ward  interstitial  diffusion  of  dissolved  molecular  water, 
which  simultaneously  reacts  with  the  Si02  network  to  form 
immobile  reaction  products.  In  order  to  accomplish  this,  we 
apply  the  technique  of  ion  implantation  to  measuring  tracer 
oxygen  diffusion  in  Si02.  Using  this  technique,  the  central 
region  within  a  thermal  oxide  layer  grown  on  silicon  is  en¬ 
riched  with  immobile  ,sO.  After  heating  in  atmospheres 
with  different  water  contents,  the  extent  of  l80  diffusion  is 
determined  by  observing  changes  in  the  concentration  pro¬ 
file  of  implanted  lsO  by  means  of  nuclear  resonance  profil¬ 
ing3  making  use  of  the  629-keV  resonance  of  the  lsO(  p,a) 
l5N  reaction.  The  observation  of  tracer  distribution  broad¬ 
ening  both  eliminates  surface  effects  and  enhances  the  sensi¬ 
tivity  of  diffusion  measurements,  enabling  tracer  diffusion 
measurements  to  be  carried  out  at  considerably  reduced 
temperatures  and  ambient  gas-phase  concentrations. 

II.  MODELS  OF  OXIDATION  AND  DIFFUSION 

Various  tracer  experiments,  cited  by  Deal  and  Grove,4 
have  established  that  the  thermal  oxidation  of  silicon  pro¬ 
ceeds  by  the  inward  transport  of  oxygen  through  the  oxide 
layer  as  a  constituent  of  molecularly  dissolved  oxygen  or 
water.  The  silicon  substrate  is  converted  to  new  oxide  by  the 
incorporation  of  atomic  oxygen  at  the  oxide-silicon  interface 
after  molecular  dissociation.  The  layer  thickens  by  forma¬ 
tion  of  new  oxide  underneath  the  existing  layer. 
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There  is  much  evidence  that  dry  thermal  oxide  layers 
formed  on  silicon  surfaces  have  the  same  structure  as  vitre¬ 
ous  silica.  X-ray  diffraction  experiments  showed  the  oxide 
layers  to  be  amorphous/'  The  density  of  the  layc  •  s  is  close  to 
that  of  fused  silica6  as  are  the  changes  in  density  induced  by 
ionizing  lcJiation.7  The  optical  properties  are  essentially 
identical,"  as  are  the  Si-O-Si  bond  angles  obtained  from  elec¬ 
tron  diffraction  and  IR  absorption  experiments.1' 

The  mechanism  of  oxygen  diffusion  through  Si02  is 
now  generally  Mieved  to  be  the  interstitial  transport  of  dis¬ 
solved  molecular  oxygen.  The  process  occurs  without  any 
direct  bonding  between  the  diffusing  species  and  the  Si02 
network;  it  rather  involves  the  transport  of  nonreacting  mol¬ 
ecules  through  the  interstices  or  cavities  in  the  network.  This 
interstitial  model  was  first  proposed  by  Anderson  and  Stu¬ 
art10;  they  calcula.  d  the  activation  energy  for  diffusion  as 
the  elastic  energy  required  to  deform  the  network  by  enlarg¬ 
ing  the  circular  doorways  between  interstices  enough  for  the 
diffusing  species  to  pass.  An  essential  feature  of  this  model  is 
a  lack  of  chemical  binding  between  the  diffusing  species  and 
the  network."  Although  the  model  is  simplistic,  the  func¬ 
tional  dependence  of  activation  energy  upon  molecular  di¬ 
ameter  agrees  well  with  activation  energies  measured  for 
molecular  diffusion  of  various  gases  in  fused  silica  (Ref.  12, 
p.  42).  It  is  also  consistent  with  the  linear  dependence  of 
diffusivity  upon  ambient  pressure  observed  by  Norton13  for 
oxygen  in  vitreous  silica,  whose  diffusion  is  not  accompanied 
by  any  significant  reactions  with  the  Si02  network.  This  was 
demonstrated  recently  in  a  tracer  experiment14  in  which 
thermal  oxide  films  which  were  grown  in  natural  oxygen 
were  further  grown  in  180-enriched  gas.  Almost  all  addition¬ 
al  ,80  was  observed  to  collect  in  a  layer  at  the  Si-SiOj  inter¬ 
face,  while  less  than  0.3%  remained  in  the  previously  grown 
region.  Great  care  was  taken  to  eliminate  any  water  in  this 
experiment. 

The  mechanism  for  water  diffusion  through  Si02  is 
similarly  thought  to  involve  the  interstitial  transport  of  dis¬ 
solved  molecular  water.2  The  transport  of  the  diffusing  spe¬ 
cies,  dissolved  molecular  water,  is  again  postulated  to  in¬ 
volve  no  direct  reaction  with  the  network.  However,  the 
permeation  of  water  in  Si02  involves  a  strong  reaction  with 
silicon-oxygen  bonds  in  the  network,  in  which  the  molecu- 
larly  dissolved  water  dissociates  to  form  two  OH  units  with 
the  additional  O  ion  being  provided  by  the  SiOj  network.2,15 
This  reversible  reaction 

HzO  +  Si-0-Si*±2  SiOH  (1) 

leads  to  equilibrium  concentrations  of  reaction  products 
which  are  large  compared  to  that  of  molecularly  dissolved 
water,  and  which  depend  on  the  square  root  of  ambient  wa¬ 
ter  vapor  concentration,  in  agreement  with  measurements  of 
water  solubility  in  vitreous  silica.16  The  reaction  products 
(i.e.,  SiOH  groups)  are  postulated  to  be  immobile;  only  the 
molecularly  dissolved  water  diffuses. 

In  the  present  study,  it  was  found  that  the  reaction  be¬ 
tween  dissolved  molecular  water  and  the  Si02  network  acts 
to  remove  previously  immobile  implanted  >sO  from  its  net¬ 
work  site  and  to  convert  it  instead  to  a  constituted  of  a 


mobile  species.  In  the  context  of  this  model,  that  species  is 
molecular  water  and  the  exchange  proceeds  via  the  reaction 

H260  +  Si-IH0-Si— >SilftOH  +  Si18OH 
— *H2"0  4-  Si-,60-Si.  (2) 

The  tracer  180  was  observed  to  diffuse  through  Si02  layers 
with  a  diffusivity  which  was  linearly  proportional  to  the  par¬ 
tial  pressure  of  ambient  water  vapor,  and  with  an  activation 
energy  of  diffusion  which  matched  those  previously  mea¬ 
sured  for  water,16  even  in  atmospheres  containing  far  larger 
concentrations  of  oxygen.  In  the  absence  of  ambient  water 
vapor,  no  diffusion  could  be  detected. 

111.  EXPERIMENTAL  PROCEDURE 

A.  Materials  preparation 

The  Si02  layers  were  formed  by  thermal  oxidation  of  2- 
in.-diam  single-crystal  silicon  wafers  obtained  from  Mon¬ 
santo  Corporation;  they  were  of  ( 100)  orientation,  with  one 
face  polished;  their  resistivity  was  4  fl  cm  {n  type).  They  were 
degreased,  cleaned,  and  oxidized  to  2000  A  in  flowing  dry 
oxygen  at  1  atm  at  1200  C.  1bO  in  the  form  of  isotopically 
enriched  gas  (70  at.  %,  obtained  from  ProChem  Isotopes) 
was  then  implanted  into  the  Si02  films  to  a  fluence  of 
3  X  1015  ions/cm2.  To  avoid  H260  contamination  the  mass 
analyzer  was  set  at  36  amu;  the  accelerating  voltage  was  80 
kV.  The  resulting  lsO  concentration  profile  was  an  Edge- 
worth  distribution  centered  at  768  A  with  a  standard  devi¬ 
ation  of  283  A  and  a  third-moment  ratio  of  —  0. 164  (Ref. 
17).  A  mass  analysis  of  the  residual  gas  within  the  ion  im¬ 
plantation  apparatus  showed  the  major  contaminants  to  be 
(160-180),  (l5N)2,  and  40 Ar,  in  quantities  insufficient  to  affect 
the  fluence  measurement  significantly.  The  wafers  received 
no  immediate  postimplant  treatment.  Samples  were  pre¬ 
pared  by  cleaving  the  oxidized  and  implanted  wafers  into 
sections  which  had  been  individually  marked  to  indicate 
their  original  locations  on  the  wafers. 

B.  Thermal  treatments 

The  thermal  treatments,  which  involved  heating  sam¬ 
ples  in  atmospheres  containing  various  water  vapor  concen¬ 
trations,  were  carried  out  in  order  to  determine  both  the 
activation  energy  of  tracer  180  diffusion  as  well  as  its  func¬ 
tional  dependence  upon  ambient  vapor  phase  concentration. 
The  former  was  determined  by  observing  the  temperature 
dependence  of  the  tracer  diffusivity  in  a  particular  atmo¬ 
sphere;  the  latter  by  observing  its  dependence  on  H20  pres¬ 
sure  at  a  single  temperature. 

1.  Temperature  dependence 

The  first  set  of  thermal  treatments  was  carried  out  in 
room  air.  The  room  was  maintained  at  a  temperature  of  2 1  C 
and  a  relative  humidity  of  50%.  Samples  were  placed  in  a 
horizontal  tubular  alloying  furnace,  lined  with  a  fused 
quartz  diffusion  tube  and  open  at  one  end;  their  tempera¬ 
tures  were  monitored  by  means  of  a  thermocouple  which 
was  inserted  in  the  diffusion  tube  at  the  sample  location. 
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Treatments  were  carried  out  by  loading  a  batch  of  samples 
on  a  preheated  fused  quartz  carrier;  sample  manipulation 
involved  withdrawing  the  carrier  from  the  furnace  for  no 
more  than  30  s  per  operation.  Temperatures  ranged  from 
413  to  789  C  and  durations  from  3540  to  230  800  >ther 
set  of  treatments  was  carried  out  in  dry  N2  at  1  atm  in  an 
identical  manner.  Gas  was  obtained  via  boiloff  of  liquid  ni¬ 
trogen;  flow  was  500  cc/min  (equivalent  at  70  ‘F,  1  atm). 
Temperatures  ranged  from  543  to  1 160  C  and  durations 
from  10  980  to  I  799  210  s.  A  similar  procedure  was  fol¬ 
lowed  for  thermal  treatments  in  steam  at  1  atm.  The  steam 
was  generated  by  boiling  deionized  water  in  a  stoppered 
flask.  The  boiling  rate  was  governed  by  a  variable  autotrans¬ 
former  supplying  power  to  a  hemispherical  heating  mantle 
around  the  flask.  Steam  was  conducted  to  the  narrowed  end 
of  the  diffusion  tube  through  a  horizontal  pyrex  tube  emerg¬ 
ing  from  the  neck  of  the  flask.  The  entire  length  of  horizontal 
tubing  and  the  neck  of  the  flask  were  wrapped  in  fiberglass 
heating  tape  and  maintained  at  a  temperature  of  120  C  to 
avoid  condensation.  An  antibackstreaming  plug  formed  by  a 
wad  of  fiberglass  tape  around  the  thermocouple  tube  al¬ 
lowed  a  slight  positive  pressure  to  be  maintained  in  the  fur¬ 
nace.  Temperatures  ranged  from  263  to  619  C  and  durations 
from  1290  to  85  400  s. 

2.  Pressure  dependence 

A  set  of  thermal  treatments  was  carried  out  in  pure 
water  vapor  at  reduced  pressures;  these  were  performed 
within  the  bell  jar  of  a  vacuum  evaporation  station.  The  sam¬ 
ples  were  placed  within  a  section  of  quartz  tubing,  as  was  a 
thermocouple.  The  tubing  was  then  inserted  into  a  small 
tubular  resistance  furnace,  whose  pedestal  mount  rested  on 
the  baseplate.  Electrical  power  for  the  furnace  was  supplied 
by  a  variable  autotransformer  through  feedthroughs  in  the 
baseplate;  roughly  80  W  was  required  to  maintain  a  sample 
temperature  of  800  C.  No  difficulty  was  encountered  in 
keeping  the  temperature  constant  to  within  ±  3  C  by  man¬ 
ually  adjusting  the  autotransformer  voltage.  Water  vapor 
was  supplied  by  evaporation  of  deionized  water  in  a  stop¬ 
pered  vacuum  flask;  a  continuous  flow  was  maintained  by 
exhausting  the  chan  ,  ir  with  a  mechanical  forepump 
through  a  choked  down  roughing  valve.  Vapor  pressure  was 
regulated  by  manually  adjusting  both  the  roughing  valve 
and  a  micrometer-type  variable  leak  valve  in  the  supply  line; 
pressure  excursions  on  the  longer  runs  amounted  to  no  more 
than  20%  of  the  average  value.  Temperatures  ranged  from 
740  to  842  C,  durations  from  7200  to  609  000  s,  and  pres¬ 
sures  from  0.062  to  8.9  Torr. 

In  addition  to  these  treatments,  several  samples  were 
treated  in  high  vacuum  ( 10~6  Torr)  in  order  to  assess  the 
substitutionality  of  the  implanted  lsO  and  the  extent  of  dam¬ 
age  to  the  oxide.  One  was  maintained  at  8 1 5  C  for  579  600  s; 
another  at  690  C  for  230  400  s,  following  which  it  received  a 
dry  nitrogen  treatment  at  1 150  C  for  1 1  400  s. 

C.  Nuclear  resonance  profiling 

In  order  to  determine  the  extent  of  180  diffusion  in¬ 
duced  by  the  thermal  treatments,  alpha  particle  yield  curves 


were  taken  using  the  method  of  nuclear  resonance  profiling.  3 
The  energy  of  a  beam  of  protons  from  a  2-MeV  Van  de 
Graaff  electrostatic  accelerator  (Model  A,  High  Voltage  En¬ 
gineering  Corp. )  was  systematically  stepped  in  the  neighbor¬ 
hood  of  the  629-keV  resonance  in  the  l80  [p,a)  l5N  reaction. 
Beam  energy  was  monitored  by  means  of  an  NMR  gauss- 
meter  whose  probe  was  mounted  between  the  analyzing 
magnet  pole  shoes;  excursions  in  the  magnetic  field  were 
easily  kept  small  enough  to  limit  beam  energy  variations  to 
under  ±  40  eV  during  a  typical  1.5-2  min  duration  alpha 
particle  count.  The  detector  was  a  silicon  surface  barrier  de¬ 
vice  (25-mm2  area,  100-//m  depletion  depth)  at  an  angle  of 
1 50°  to  the  beam  direction,  shielded  from  backscattered  pro¬ 
tons  by  a  1.3-mg  cm  2  polyester  foil.  Samples  were  mounted 
normal  to  the  beam  and  biased  to  collect  secondaries.  Each 
batch  of  samples  profiled  included  a  reference  sample  that 
had  been  lsO  implanted  but  not  thermally  treated;  the  refer¬ 
ence  samples  provided  a  continuous  self-calibration  for  the 
proton  beam  energy.  An  additional  calibration  was  per¬ 
formed  by  means  of  ( p,  y)  nuclear  resonance  profiling  of 
thick  targets  containing  uniform  concentrations  of  appropri¬ 
ate  nuclides.  This  was  done  primarily  to  determine  the 
spread  (i.e.,  the  distribution  of  proton  energies)  in  the  inci¬ 
dent  beam. 

The  alpha  particle  counts  at  each  beam  energy  accumu¬ 
lated  on  a  scaler  until  the  proton  fluence  on  the  sample  had 
reached  a  preset  value  (typically  50 pC).  Alpha  particle 
counts  ranged  from  1000  to  5000  per  run;  background 
counts  were  negligible.  Beam  current  was  held  to  0.5  pA  to 
avoid  heating  the  samples.  The  range  of  beam  energy  settings 
spanned  was  typically  1 3  keV  on  either  side  of  the  energy  at 
which  the  alpha  yield  was  maximum,  corresponding  to  a 
distance  of  2 1 00  A  on  either  side  of  the  distribution  peak. 1 8 

IV.  RESULTS  AND  ANALYSIS 

The  raw  data  obtained  by  the  nuclear  resonance  profil¬ 
ing  consisted  of  an  alpha  particle  yield  curve  for  each  sam¬ 
ple.  A  representative  set  of  these  curves  appears  in  Fig.  1,  in 
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FIG.  1.  Observed  alpha 
particle  yield  curves  for  se¬ 
lected  samples.  Continuous 
curves  represent  calculated 
yield  curves. 
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which  the  different  symbols  denote  data  recorded  for  the 
samples  indicated.  For  the  two  lowest  curves  the  counting 
time  was  doubled  to  obtain  sufficient  statistical  accuracy:  the 
points  plotted  for  those  curves  are  half  the  recorded  counts. 
Also  shown  with  each  observed  curve  is  a  theoretically  de¬ 
rived  curve  (the  continuous  curve  between  the  data  points); 
these  were  calculated  according  to  the  procedure  described 
in  the  following  section. 

A.  Method  of  determining  iaO  tracer  dlffuslvity 

As  discussed  below,  the  implantation  process  caused 
l80  to  effectively  replace  existing  oxygen  as  a  constituent  of 
the  Si02  network.  According  to  the  model  of  diffusion  de¬ 
scribed  in  Sec.  II,  the  exchange  of  network  oxygen  with  dis¬ 
solved  molecular  water  is  solely  responsible  for  any  observed 
change  in  the  distribution  of  implanted  1  sO.  The  dependence 
of  the  tracer  ,80  diffusion  on  thermal  treatment  parameters 
is  a  diagnostic  of  the  mechanism  of  water  diffusion  in  Si02 
layers,  and  a  test  of  the  validity  of  this  model. 

The  analysis  of  tracer  ,sO  transport  via  the  exchange 
reaction  involves  three  interrelated  diffusion  coefficients. 
The  first  of  these,  Dt„,  is  the  quantity  normally  measured  in 
permeation  experiments.  This  effective  diffusivity  is  inde¬ 
pendent  of  ambient  gas-phase  concentration  only  if  the 
fluxes  of  emerging  particles  and  dissolved  species  have  the 
same  pressure  dependence.  (Such  is  not  the  case  for  water 
diffusion,  where  Dcir  is  the  apparent  diffusivity  of  reacted 
water.)  The  transport  of  the  actual  mobile  species  is  de¬ 
scribed  by  the  microscopic  diffusivity  D.  For  a  dilute  gas  in  a 
spatially  uniform  network,  D  is  independent  of  ambient  gas- 
phase  concentration,  spatial  location,  and  dissolved  gas  con¬ 
centration;  it  depends  only  on  temperature.  For  oxygen  dif¬ 
fusion  in  Si02,  which  proceeds  without  appreciable  dissocia-. 
tion  or  reaction  with  the  network,  Dcb  and  D  have  the  same 
value.  However,  for  water  diffusion  they  are  related  by  the 
equation 

D"  =Z)2[H2O]/(2[H20]  +  [  SiOH] ).  (3) 

The  fraction  represents  the  probability  that  any  particular 
OH  group  will  be  part  of  a  diffusing  water  molecule;  the 
factor  of  2  accounts  for  the  removal  of  two  OH  groups  from 
the  network  by  each  water  molecule.  The  effective  diffusion 
of  tracer  ,80  by  means  of  the  exchange  reaction  is  described 
by  the  tracer  diffusivity  D  *.  In  the  presence  of  free  exchange 
under  quasisteady  conditions,  D  *  is  related  to  D  by  the 
equation 

D*  =  D  [H20]/(2[Si02]  +  [H20]  +  [SiOH]).  (4) 

The  fraction  is  the  probability  that  an  ,sO  will  be  part  of  a 
diffusing  water  molecule.  The  factor  of  2  accounts  for  the 
incorporation  of  two  oxygen  atoms  per  Si02  molecule. 

The  value  of  D  *  characteristic  of  a  given  thermal  treat¬ 
ment  is  reflected  in  the  change  of  the  implanted  l80  distribu¬ 
tion  resulting  from  the  treatment.  Once  the  alpha  particle 
yield  is  measured,  the  extent  of  change  is  in  principle  deter¬ 
minable  by  inversion  of  the  integral  equation  for  the  yield, 
namely, 


Y(E„)=  ['dxpix^r  dEig{Eb^) 

Jo  Jo 

X  f  ‘dE^-(E)f{E,Eis),  (5) 

Jo  all 


where  Eb  is  the  mean  energy  of  the  incident  proton  beam,  E 
is  the  energy  of  reacting  protons,  E,  is  the  normal  energy  of 
protons  incident  on  the  surface  of  the  layer,  x0  is  the  thick¬ 
ness  of  the  oxide  layer,  x  is  the  depth  at  which  the  reaction 
occurs,  p(x,t )  is  the  lsO  concentration  profile  (i.e.,  atoms  per 
unit  volume  at  depth  x  at  time  t ),  do/dfl  is  the  differential 
cross  section  of  the  ,80  (p,a)  l5N  reaction,  g{Eh  ,El  )dE,  is  the 
fraction  of  protons  in  the  beam  which  have  an  incident  be¬ 
tween  Ej  and  E,  -f  dE,,  and  f[E,E,pc)dE  is  the  probability 
that  a  proton  which  has  started  with  energy  E,  and  arrived  at 
depth  x  will  have  an  energy  between  E  and  E  +  dE. 

If  the  functions  g(Eb  ,E,  ),/(£,£,  ,*),  and  da/dO  ( E )  are 
all  previously  known  and  the  yield  Y (Eb)  is  observed,  then 
the  function p(x)  can  in  principle  be  determined  by  inversion 
of  the  integral  equation.  An  enlightening  discussion  of  this 
matter  appears  in  Ref.  19,  in  which  a  similar  calculation 
[involving  the 27 Al(  p,y)  28Si  reaction]  was  done. 

In  practice,  deconvolution  is  a  complicated  procedure, 
entailing  large  errors  whose  propagation  is  difficult  to  track. 
However,  once  the  initial  l80  distribution  is  determined,  its 
subsequent  evolution  is  predictable  through  macroscopic 
diffusion  theory  provided  the  properties  of  the  medium  are 
correctly  perceived  and  the  proper  boundary  conditions  are 
applied.  This  is  treated  further  in  the  Appendix,  where  the 
evolved  distribution  is  calculated  as  a  functional  of  the  prod¬ 
uct  D  *t  [t  being  the  duration  of  treatment).  Based  on  this 
distribution  together  with  the  other  constituent  functions  of 
the  yield  integral  given  in  the  Appendix,20  the  yield  was  cal¬ 
culated  for  each  of  a  number  of  discrete  values  of  D  *t.  The 
resulting  family  of  yield  curves,  which  is  illustrated  in  Fig.  2, 
was  compared  to  the  experimental  yield  curve  resulting  from 
that  thermal  treatment.  The  appropriate  value  of  D  *t  was 


FIG.  2.  Calculated  alpha  particle  yield  curves.  Labels  indicate  D  */,  in  units 
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taken  as  that  which  indexed  the  calculated  yield  curve  most 
closely  matching  the  experimental  yield  curve. 

For  routine  measurements  of  D  *,  a  procedure  was 
adopted  which  allowed  shorter  thermal  treatments  and  also 
greatly  shortened  the  time  spent  in  nuclear  profiling.  It  can 
be  seen  from  the  family  of  theoretical  yield  curves  in  Fig.  2 
that  the  peak  yield  of  alpha  particles  can  be  regarded  as  a 
monotonically  decreasing  function  of  D*t  This  function  was 
constructed  by  calculating  a  large  number  of  theoretical  al¬ 
pha  particle  yield  curves.  The  peak  value  of  each  of  those 
curves  was  tabulated  along  with  the  vaiue  of  D  *t  associated 
with  the  curve;  those  values  were  then  connected  with  a 
smooth  curve.  To  determine  the  proper  value  for  a  given 
thermal  using  this  function,  it  was  necessary  to  measure  an 
alpha  particle  yield  curve  for  a  sample  only  in  the  vicinity  of 
its  peak  so  as  to  determine  the  peak  yield.  This  quantity  was 
then  normalized  to  that  of  a  reference  sample  located  nearby 
on  the  wafer.  The  appropriate  value  of  D  *t  was  then  found 
simply  by  locating  the  resulting  peak  yield  ratio  on  a  graph  of 
that  function. 


B.  Temperature  and  pressure  dependence  of  tracer 
diffusivity 

The  temperature  dependence  of  D*  was  determined  as 
follows.  First,  D  *  was  calculated  from  the  observed  data  us¬ 
ing  the  methods  described  in  the  previous  section.  For  each 
of  the  three  atmospheres  (i.e.,  steam  at  1  atm,  room  air,  and 
flowing  dry  nitrogen)  the  values  of  D  *  were  fit  by  an  equation 
of  the  Arrhenius  type  with  a  temperature-independent 
preexponential  factor.  The  values  of  D  *  were  taken  as  mea¬ 
sured:  they  wer ;  not  adjusted  for  the  temperature  variation 
of  gas  phase  concentration.  The  fit  was  carried  out  by  the 
method  of  unweighted  linear  least  squares  using  \/T  and 
log ,,/)  *  as  the  independent  and  dependent  variables,  respec¬ 
tively.  The  results,  which  appear  in  Fig.  3  along  with  the 
individual  values  of  D  *,  are 

D  *Mm  =  2.7  X10  10  ±  0Me  ~ 116  9  ±  13)kc*'/rn°'/*7-crn7s,  (6) 
D*^  .ir  =  1.3xlO-12±  0  34e-|l55±14|k“l/n,ol/',rcmVs,  (7) 

/>*,  N ,  =  1.8X10  'V  ,V7kc.l/mol/*7cm2/s  (g) 


TEMPERATURE  - C 


FIG.  3.  Observed  temperature  dependence  of'*0  tracer  diffusivity  at  I  atm. 
Solid  lines  represent  linear  least -squares  fits  to  the  data  shown. 
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FIG.  4.  Observed  pressure  dependence  of  '*0  tracer  diffusivity  in  water 
vapor  at  820  C.  Labels  indicate  pressure  in  fractions  of  1  atm. 

The  indicated  errors  are  derived  solely  from  the  variances  of 
the  fitting  parameters.  The  standard  deviations  in  the  fits  for 
log/)  *  are  0.232  and  0. 1 35  for  steam  and  room  air,  respec¬ 
tively,  which  imply  relative  errors  of  171%  and  136%  for 
D  * .  The  dry  nitrogen  data  points  fell  fortuitously  close;  their 
error  was  not  calculated.  The  three  fitted  lines  shown  in  Fig. 
3  are  almost  parallel;  at  600  C  that  representing  D  *Mtn  ,lr  lies 
a  factor  of  93  below  D  ,  while  D  N>  lies  a  factor  of  8. 1 
below  that. 

The  pressure  dependence  of  D  *  was  determined  by  cal¬ 
culating  D  *  in  a  similar  manner  for  the  samples  heated  at  low 
pressure.  As  they  had  been  heated  at  slightly  different  tem¬ 
peratures  their  values  of  D  *  were  adjusted  to  a  common  tem¬ 
perature  (of  820  C)  by  assuming  a  temperature  dependence 
of  that  seen  for  steam.  The  resulting  values  of  D  *  (per  unit 
pressure)  are  displayed  in  Fig.  4,  as  is  D  at  1  atm  which 
was  calculated  using  Eq.  (6).  The  indicated  errors  are  those 
arising  purely  from  counting  statistics.  As  can  be  seen,  these 
values  are  consistent  with  a  linear  pressure  dependence  over 
the  entire  range  of  pressures  studied. 

In  addition  to  these,  measurements  were  made  on  the 
samples  which  had  undergone  high-vacuum  treatments  in 
order  to  assess  substitutionality  of  the  implanted  l80  and  its 
damage  to  the  oxide.  No  difference  could  be  discerned  be¬ 
tween  the  alpha  particle  yield  obtained  from  the  sample 
which  had  been  maintained  at  815  C  for  579  600  s  at  10~6 
Ton  and  that  from  an  untreated  reference  sample.  Another 
sample  had  its  alpha  particle  yield  measured  twice.  After 
being  maintained  at  690  C  for  230  400  s  at  10~6  Torr,  the 
yield  was  compared  to  that  of  an  untreated  reference  sample; 
again  no  difference  could  be  detected.  It  was  then  main¬ 
tained  for  10  980  s  at  1 1 50  C  in  flowing  dry  nitrogen  at  1  atm 
and  the  yield  was  again  compared  to  that  of  reference  sam¬ 
ple;  this  time  the  results  implied  a  value  of  D  *  close  to  that  of 
a  sample  which  had  undergone  a  similar  treatment  in  dry 
nitrogen  but  had  no  prior  vacuum  treatment. 

V.  DISCUSSION 

A.  Absence  of  Ion  bombardment  Induced  damage  In 
StO,  network 

To  be  able  to  compare  the  diffusivity  measurements 
performed  in  this  study  to  others  made  in  Si02  films  and 
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membranes,  it  was  necessary  to  verify  that  the  properties  of 
the  material  used  here  match  those  of  the  others.  For  this 
purpose,  evidence  was  collected  which  indicated  that  the  ,80 
implantation  does  not  produce  enough  network  damage  to 
appreciably  alter  the  oxide’s  structure  or  properties,  at  least 
as  far  as  l80  tracer  diffusion  is  concerned. 

First,  the  amount  of  l80  implanted  was  small  enough  to 
avoid  gross  morphological  changes  brought  about  by  its  in¬ 
corporation,  such  as  blister  or  bubble  formation.  (None  were 
seen.)  The  actual  quantity  of  additional  oxygen  implanted  in 
the  network  amounted  to  at  most  4.2  X  1020  molecules/cm3, 
which  is  1.9%  of  the  concentration  of  02  in  Si02.21 

Second,  the  shape  of  every  observed  alpha  particle  yield 
curve  could  be  well  matched  to  one  or  another  of  the  set  of 
theoretical  curves.  These  were  calculated  on  the  basis  of  a 
dilute  homogeneous  medium;  there  was  therefore  no  indica¬ 
tion  ever  observed  of  a  region  of  anomalous  diffusivity. 

Third,  similar  thermal  treatments  (in  dry  N2  at  1 1  SO  C) 
were  performed  on  two  samples,  one  of  which  had  been 
vacuum  annealed  at  690  C  for  230  400  s.  The  l80  tracer 
diffusivity  was  found  to  be  totally  unaffected  by  this  anneal. 
These  observations  can  be  interpreted  in  the  context  of  pre¬ 
vious  investigations  of  annealing  of  radiation  damage  in 
Si02. 22-25  These  studies  indicate  that  whatever  damage  is 
produced  in  SiOz  is  entirely  removed  by  a  900  C  inert  atmo¬ 
sphere  or  vacuum  anneal.26  It  is  reasonable  to  conclude  that 
the  thermal  treatments  performed  here  to  induce  lsO  diffu¬ 
sion  were  by  themselves  sufficient  to  anneal  out  any  damage 
induced  by  the  implantation. 

To  verify  that  the  proton  bombardment  involved  in  the 
nuclear  resonance  profiling  did  not  induce  any  alteration  of 
the  network  or  diffusion  of  180,  alpha  particle  yield  curves 
were  collected  in  untreated  reference  samples  both  before 
and  after  lengthy  proton  bombardments  (about  2000 /zC, at  a 
current  of  0.5  /zA).  No  difference  could  be  detected. 

Finally,  to  determine  whether  the  implanted  180  actu¬ 
ally  replaced  existing  oxygen  in  the  network,  or  remained 
interstitial  and  free  to  migrate,  an  alpha  particle  yield  curve 
was  taken  for  an  untreated  reference  sample  after  a  high- 
vacuum  anneal  at  815  C  for  579  600  s.  If  an  appreciable 
amount  of  l80  had  been  mobile,  then  it  would  have  diffused 
during  that  period  with  a  diffusivity  characteristic  of  oxygen 
at  that  temperature,  namely,  10-9  cm2/s  (Ref.  13).  A  calcu¬ 
lation  showed  that  the  thermal  treatment  given  the  sample 
was  enough  to  guarantee  the  escape  of  virtually  all  such  mo¬ 
bile  ,80  through  the  oxide  surface.  However,  the  observed 
yield  curve  was  indistinguishable  from  that  of  an  untreated 
unannealed  reference  sample,  indicating  no  shift  in  the  ini¬ 
tial  l80  concentration  profile.  Both  curves  closely  matched 
the  theoretically  predicted  curve  for  such  an  implant,  dem¬ 
onstrating  that  the  initial  distribution  of  180  formed  by  the 
implantation  process  consists  of  bound  network  oxygen  with 
an  insignificant  fraction  remaining  mobile.  This  conclusion 
is  in  accord  with  studies  of  A1  implant  substitutionality  in 
Si02  (Ref.  27),  as  well  as  with  Primak’s  thermal  spike  model 
of  heavy  ion  energy  loss,28  which  predicts  bond  cleavage  by 
thermal  spikes  involving  peak  temperatures  of  1 1 50  C  affect¬ 
ing  3 1  000  network  oxygen  atoms  per  ion,  on  a  time  scale  of 
10~ 12  s;  rapid  refreezing  ensures  almost  complete  incorpora¬ 


tion  of  implanted  l80  as  a  constituent  of  the  network 


B.  1(0  tracer  diffusion 

Once  ambient  water  vapor  was  introduced,  the  im¬ 
planted  ,80  (which  was  found  to  be  immobile  in  vacuum  at 
temperatures  exceeding  800  C)  was  seen  to  move  in  a  manner 
consistent  with  its  exchange  with  dissolved  molecular  water 
and  its  subsequent  diffusion  as  a  constituent  of  that  species. 
This  conclusion  is  supported  by  the  following  observations: 

(i)  The  observed  linear  dependence  of  the  l80  tracer  dif¬ 
fusivity  upon  ambient  vapor  phase  concentration  is  identical 
to  that  expected  for  dissolved  molecular  water.  This  linear 
dependence  extended  over  four  orders  of  magnitude  (see  Fig. 
4),  from  1  atm  to  79  ppm  (of  1  atm),  in  pure  water  vapor — no 
carrier  gases  were  employed.  When  room  air  was  used  as  a 
carrier  gas,  the  observed  tracer  diffusivity  conformed  to  this 
linear  dependence,  but  upon  the  partial  pressure  of  atmo¬ 
spheric  water.  This  can  be  seen  from  Fig.  3,  where  the  mea¬ 
sured  ratio 


D*  /D*  —  93 

is  close  to  the  ratio  of  water  vapor  partial  pressures  for  these 
atmospheres,  namely  87  (assuming  room  air  at  20  C  and 
50%  relative  humidity).  And  when  samples  were  treated  in 
high  vacuum,  no  l80  diffusion  could  be  detected  even  after 
6.7  days  at  815  C. 

(ii)  The  temperature  dependence  of  D  *„m  showed  an 
activation  energy  of  16.9  ±  1.3  kcal/mol,  which  is  indistin¬ 
guishable  from  those  previously  observed  for  water  diffusion 
either  by  permeation  in  vitreous  silica16  or  by  growth  rate 
kinetics  of  thermal  oxidation.4-29  These  range  from  16  to  18 
kcal/mol.  This  value  is  clearly  inconsistent  with  those  pre¬ 
viously  observed  for  oxygen  diffusion, 4-,3',4-29-3°  the  most  re¬ 
liable  of  which  are  about  28  kcal/mol. 

(iii)  The  magnitude  of  D  *  is  close  to  that  expected  for 
free  exchange  between  network  oxygen  and  dissolved  water. 
In  the  presence  of  free  exchange,  D  *  is  related  to  D",  the 
experimentally  measured  diffusivity  of  water,  according  to 
Eqs.  (3)  and  (4)  by 

{D*/Dca)fmexchmtt  =  15X10”4  at  600  C 
=  7.5X  10-4  at  1200  C. 


This  ratio,  whose  value  is  derived  from  the  water  solubility 
measurements  of  Moulson  and  Roberts  (M-R),16  can  also  be 
evaluated  by  using  the  results  for  D*c,m  obtained  in  this 
study.  According  to  Eq.  (6),  this  is 


flit,., (this  work) 
D"(  M-R) 


=  6.0X  10“4  at  600  C 
=  4.4X  10~4  at  1200  C. 


The  magnitude  of  D  * _ is  thus  seen  to  be  consistent  with 

free  exchange  at  both  temperatures  measured  by  Moulson 
and  Roberts  to  within  a  factor  of  2-3,  which  is  within  the 
limits  of  accuracy  obtained  for  the  D  *  measurements  in  this 
study. 

(iv)  In  the  presence  of  atmospheric  oxygen,  the  activa¬ 
tion  energy  for  ,jr  remained  that  of  water  [roughly  16 
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kcal/mol,  as  indicated  in  Eq.  (7)],  rather  than  reverting  to 
that  of  oxygen  (28  kcal/mol).  This  was  true  even  at  800  C, 
where  D(h  exceeds  (£>„«- )H,o  by  an  order  of  magnitude,  and 
even  though  the  concentration  of  atmospheric  oxygen  was 
1 7  times  that  of  water  vapor  (20%  of  760  Torr  vs.  50%  of  1 8 
Torr)  This  implied  the  absence  of  exchange  between  atmo¬ 
spheric  and  network  oxygen,  in  agreement  with  similar  re¬ 
cent  findings  by  Rosencher  et  al ., M  who  observed  the  failure 
of  diffusing  l802  to  replace  existing  network  l60  during  ther¬ 
mal  growth  of  oxide  films  on  silicon  in  dry  lsO  atmospheres. 

In  addition  to  supporting  the  above  model  of  water  dif¬ 
fusion  in  Si02,  the  measurements  made  in  this  study  pro¬ 
vided  evidence  indicating  that  chemical  equilibrium  be¬ 
tween  water  dissolved  in  Si02  films  and  its  reaction  products 
occasionally  had  not  been  attained  over  the  time  spans  of  the 
thermal  treatments.  The  values  of  D  *  scattered  in  a  seeming¬ 
ly  systematic  manner:  the  smallest  values  tended  to  be  asso¬ 
ciated  with  the  shortest  thermal  treatments  as  well  as  with 
the  performance  of  a  high-vacuum  anneal  prior  to  thermal 
treatment,  which  was  done  for  some  samples.  Further,  the 
amount  of  yield  curve  broadening  in  these  cases  was  anoma¬ 
lously  small  compared  to  the  decay  in  peak  yield.  This  same 
phenomenon  seems  to  have  been  observed  in  some  early 
studies  of  l80  permeation  in  Si02  fibers  and  membrances.31 
A  more  thorough  study  of  the  effects  of  nonequilibrium  per¬ 
meation  conditions  on  the  evolution  of  ,80  distributions  is 
currently  in  progress. 

VI.  CONCLUSIONS 

The  results  of  the  lsO  tracer  diffusion  measurements 
performed  in  this  study  clearly  demonstrate  that  water  diffu¬ 
sion  in  Si02  films  proceeds  by  the  transport  of  molecular 
water  and  that  it  is  accompanied  by  a  strong  reversible  reac¬ 
tion  between  molecular  water  and  network  oxygen.  The  ob¬ 
served  linear  dependence  of  tracer  diffusivity  upon  partial 
pressure  of  atmospheric  water  vapor,  for  pressures  ranging 
from  1  atm  to  as  little  as  0.06  Ton,  indicates  that  network 
oxygen  is  exchanged  with  molecular  water.  The  activation 
energy  of  tracer  diffusion  remained  about  1 6  kcal/mol,  char¬ 
acteristic  of  water  diffusion,  over  a  range  of  temperatures 
from  1  I  f  0  to  260  C.  Further,  the  magnitude  of  the  tracer 
diffusivity  was  consistent  with  the  presence  of  free  exchange 
between  network  oxygen  and  diffusing  water.  The  extent  of 
tracer  diffusion  proved  to  be  unaffected  by  the  presence  of 
atmospheric  oxygen,  indicating  the  absence  of  exchange 
with  that  species.  The  results  obtained  here  agreed  with  cor¬ 
responding  results  obtained  by  most  other  investigators  us¬ 
ing  different  means.  They  are  in  complete  agreement  with 
recently  proposed  models  of  oxygen  and  water  diffusion  in 
Si02  in  which  the  diffusion  proceeds  by  transport  of  the  mo¬ 
lecular  species  through  interstices  in  the  network  without 
any  direct  reaction  between  the  molecular  species  and  the 
Si02  (the  strong  exchange  reaction  accompanying  water  dif¬ 
fusion  involves  immobile  OH  groups  which  are  dissociation 
products  of  the  molecular  water). 

The  process  employed  for  introducing  tracer  180  into 
the  network,  that  of  ion  implantation,  was  shown  to  produce 
a  region  within  the  oxide  in  which  the  180  was  fully  incorpo- 
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rated  in  the  network:  no  l80  diffusion  was  observed  in  high 
vacuum,  even  after  161  h  at  815  C.  The  properties  of  the 
network  remained  otherwise  unaltered,  at  least  as  far  as  dif¬ 
fusion  was  concerned.  This  was  demonstrated  both  by  the 
close  match  of  the  data  to  theoretical  predictions  based  on  a 
dilute  homogeneous  medium  as  well  as  by  the  absence  of  any 
effect  on  diffusion  of  a  high-temperature  treatment  which 
was  sufficient  to  anneal  out  any  implantation-induced  dam¬ 
age.  These  findings  indicate  that  the  thermal  treatments  per¬ 
formed  to  induce  180  diffusion  were  themselves  sufficient  to 
anneal  out  any  damage. 

The  technique  of  nuclear  resonance  profiling  proved  to 
be  a  convenient,  nondestructive,  and  sensitive  method  of 
measuring  l80  diffusivities  as  low  as  10“ 17  cmVs  at  tem¬ 
peratures  as  low  as  260  C. 
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APPENDIX.  CONSTITUENT  FUNCTIONS  OF  THE  ALPHA 
PARTICLE  YIELD  INTEGRAL 

p(x,t ):  The  Si02  film  is  assumed  to  be  a  dilute  homogen¬ 
eous  medium  with  respect  to  water  diffusion;  the  concentra¬ 
tions  of  dissolved  molecular  water  and  reacted  water  remain 
at  their  equilibrium  levels.  If  the  film  is  heated  for  a  duration 
l  under  those  conditions,  then  the  distribution p[x,t )  of  180  in 
excess  of  the  natural  background  level  of  0.02%  isotopic 
fraction  will  change  with  time  so  as  to  satisfy  the  diffusion 
equation 

p{x,t)  =  Dt~p{x,t), 
ax 

with  D  *  independent  oip,  x,  and  t.  The  lsO  is  assumed  to  be 
transported  only  as  a  constituent  of  molecularly  dissolved 
water.  The  boundary  conditions  are  therefore  determined  by 
considering  the  behavior  of  this  substance  at  the  boundaries. 

(i)  At  the  oxide  surface,  equilibrium  between  the  outside 
atmosphere  and  water  dissolved  in  the  oxide  occurs  rapidly; 
surface  absorption  and  reactions  are  not  involved  (Ref.  1 2,  p. 
128).  This  implies  the  condition  p[0,t )  =  0. 

(ii)  At  the  oxide-silicon  interface,  the  boundary  condi¬ 
tion  is  temperature  dependent.  This  can  be  seen  from  the 
analysis  of  steady-state  steam  oxidation  performed  by  Deal 
and  Grove.4  According  to  their  Eqs.  (6)  and  (7),  the  ratio  of 
molecularly  dissolved  water  concentration  at  the  interface  to 
that  at  the  outer  surface  is  given  by 

Si.  =  i 

Co  1  +  kXo/D^g 

where  k,  the  interface  rate  constant,  is  determined  through 
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observations  of  the  linear  and  parabolic  rate  constants  .4  and 
B  and  the  equilibrium  solubility  C  *  (x0  is  the  layer  thickness). 
Deal  and  Grove  give  a  value  of 

k  =  1.8X  103  /itn/h 

at  1000  C  with  an  activation  energy  of  45.3  kcal/mol.  Using 
this  in  conjunction  with  the  Moulson  and  Roberts 
equation16 

D"  =  10~6exp[  -  1 8.3(kcal/mol)/R T ] ,  (cmVsec) 

the  temperature  dependence  of  C,/C0  is  seen  to  be  the 
following: 


T(C) 

kXo/D't 

Ct/C0 

250 

3.1X10-7 

~1 

475 

7.7X10”4 

~1 

600 

0.011 

0.989 

1000 

1.4 

0.417 

1200 

5.9 

0.145 

At  temperatures  encountered  in  this  work,  the  ratio  C,/C0  is 
more  than  98%.  The  appropriate  boundary  condition  is 
therefore  dp(x0,t  )/dx  =  0. 

(iii)  For  sufficiently  long  times  all  the  excess  lsO  will 
eventually  disappear,  either  through  surface  leakage  or  in¬ 
terface  incorporation  of  diffusing  molecular  water,  implying 
the  condition  p(x,  oo )  =  0. 

(iv)  The  initial  distribution  p(x, 0)  was  taken  to  be  that 
implantation  profile  given  by  Gibbons  et  al.17  for  the  appro¬ 
priate  experimental  conditions.  This  was  an  Edgeworth  dis¬ 
tribution  with  a  mean  projected  range  of  768  A,  a  standard 
deviation  of  283  A,  and  a  third  moment  ratio  of  —  0. 163. 
The  solution  to  the  diffusion  equation  is  straightforward, 
being 

p[x,t )  =  V  B2n_te~  '<2" -  'K^.)i1"*'sin  ( ^ x] , 

n  =  I  V  2  X0  / 

with  , 

Bk  =  2  I  dx  sin(J/cjrx^(jc,0). 

Jo 

g[EbJZi):  The  distribution  of  energies  in  the  incident 
beam  was  found  experimentally  through  analysis  of  the 
thick-target 27 Al (p,y)  28Si  excitation  curve.  The  energy 
spread  amounted  to  less  than  1  keV;  the  shape  was  taken  to 
be  a  Gaussian  with  a  standard  deviation  of  1  keV. 

f(E,EjPc):  After  some  test  calculations  with  Vavilov  dis¬ 
tributions32  it  was  found  that  for  the  particular  conditions  of 
this  experiment  a  sufficiently  accurate  approximation  to  the 
energy  straggling  distribution  was 

/M  )  =  (2 irk,x)  ~  1/2exp[  -  (A  -  e x)2/2ksx], 

where  A  =  Et  —  E(x),  k,  =  17.75  keV2/im_l  (Ref.  33),  and 
e  =  61  keV  /jm-1  (Ref.  18). 

da/dfl  [E ):  The  recent  data  of  Amsel  et  al.  (Ref.  34,  p. 
163)  in  the  region  of  the  629-keV  resonance  was  well  fit  by  a 
sum  of  two  interfering  Lorentzians  plus  an  exponential 
background: 

da  \r2m(\2Ji-K)-{E-Em  )4.7e|g  ~ 

dfl  {E-E^  +  mj 

+  Ke[E  ~  £r”|/!8mb/steradian, 


where 

K  =  2  ( E<ErJ 

=  1.75  ( E>ErJ , 
rm  =  1.9  keV, 

E,a  =  629  keV. 
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